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Abstract 

Numerical simulations of the current-voltage characteristics of an ordered 
two-dimensional Josephson junction array at an irrational flux quantum per 
plaquette are presented. The results are consistent with an scaling analy- 
sis which assumes a zero temperature vortex glass transition. The thermal 
correlation length exponent characterizing this transition is found to be sig- 
nificantly different from the corresponding value for vortex-glass models in 
disordered two-dimensional superconductors. This leads to a current scale 
where nonlinearities appear in the current-voltage characteristics decreasing 
with temperature T roughly as T'^ in contrast with the T^ behavior expected 

for disordered models. 
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There has been an increasing interest in two-dimensional Josephson-j unction arrays as 
a model system to study disorder and frustration effects as found in high-Tc superconduc- 
tors and spin-glass systems^. In Josephson-j unction arrays, frustration without disorder 
can in principle be introduced by applying an external magnetic field on a perfect periodic 
array. The frustration parameter /, the number of flux quantum per plaquette, is given 
by / = 0/0O, the ratio of the magnetic flux through a plaquette to the superconducting 
flux quantum 0o = hc/2e, and can be tuned by varying the strength of the external field. 
Frustration can be viewed as resulting from a competition between the underlying periodic 
pinning potential of the array and the periodicity of the vortex latticei. At a rational value 
of /, the ground state is a commensurate pinned vortex lattice leading to discrete symme- 
tries in addition to the continuous U{1) symmetry of the superconducting order parameter. 
In particular, for / = 1/2 which has been intensively studied both experimentally and the- 
oretically, a superconducting phase transition takes place at finite temperatures with an 
interplay of f/(l) and discrete Z2 symmetryi. At irrational values of /, the behavior is much 
less understood since the vortex lattice is now incommensurate with the array. The ground 
state consists of a disordered vortex pattern^ lacking long range order which can also be re- 
garded as a vortex-glass state without disorder. One can not completely rule out a possible 
glass transition at finite temperatures as has been suggested by Halseyi but other arguments 
suggest a zero temperature transitionii. In any case, on the basis of a close analogy between 
this system and gauge glass models of disordered superconductors, one expects metastable 
states with long relaxation times and a nonzero critical current at zero temperatureBi. In 
fact, measurements of I-V characteristics in two-dimensional superconducting wire networks 
at an irrational frustrationi have been interpreted within the scaling analysis of a vortex 
glass transition as in three-dimensional disordered superconductors where there is evidence 
of a finite-temperature glass transition! In the case of superconducting wires, the behavior 
is likely to be dominated by a mean field transition which should correspond to the behavior 
at higher dimensions. This could provide a possible explanation for the observed scaling 
behavior of a finite-temperature transition. On the other hand, no evidence was found in 



experiments on two-dimensional proximity-coupled Josephson-j unction arrays! where phase 
fluctuations are expected to be more important. In disordered superconductors, studies of 
different models of the vortex glass tend to agree that in two dimensions a vortex glass transi- 
tion takes place only at zero temperatur J which is supported both by numerical simulation J 
and experiment^^. Although being completely different in the nature of their ground states, 
the similarity of the behavior of these two systems, specially regarding slow relaxation vortex 
dynamics and a possible zero temperature transition in low dimensions, strongly suggests 
that a two-dimensional array at irrational / should behave as a zero-temperature vortex 
glass. In this case, the appropriate scaling of the I-V characteristics should be the one cor- 
responding to a zero temperature transition while the different nature of the ground states 
should be reflected in the value of the critical exponents as a different universality class. It 
seems therefore worthwhile to investigate to which extent an array at irrational / can be 
described as a zero temperature vortex glass. 

In this work, we present simulations of the current-voltage characteristics of a Josephson 
junction array at an irrational flux quantum per plaquette / = (3 — v^)/2, a golden irra- 
tional, and an scaling analysiJ^ which assumes a zero temperature vortex glass transition. 
The results are consistent with the scaling assumption and allow for an estimation of the 
thermal correlation length critical exponent v characterizing the zero temperature transi- 
tion. This critical exponent is found to be signiflcantly different from the corresponding 
value for vortex-glass models in disordered two-dimensional superconductors. As a result, 
the current density scale, Jni ~ T^^*^, where nonlinearities appear in the current-voltage 
characteristics decreases with temperature roughly as in contrast with the behavior 
expected for disordered mo deJ. This could provide a signature of the behavior at irrational 
frustration in experimental studies of ordered arrays of Josephson junctions. 

We considered an ordered array of Josephson junctions deflned on a square lattice. The 
I-V characteristics of the array was computed using an overdamped Langevin molecular- 
dynamics as described by Falo et al0. One allows for a capacitance to the ground Co in 
addition to a shunt resistance i?o between superconducting grains. In the overdamped hmit, 



this particular dynamics reduces to the standard resistance shunt junction (RSJ) modelEl 
commonly used in dynamical simulationslll. The Langevin equations can be written as 

+ ^ E ^-^^^ = E M0^ - 0, - A.) + ir + E (1) 
° j j j 

where rjij is a Gaussian white noise satisfying 

< r]ij{t) >= 0, < Vijit)Vkiit') >= —^6ij^kiS{t - t') (2) 

and j are nearest-neighbor pairs. is the external current and Ic the Josephson-j unction 
critical current. The bond variables Aij are constrained to Y^ij = 27r/ where represent 
a direct sum around a plaquette. We use units where h/2e = l,Ro = 1, Ic = 1 and set the 
parameter I^R^^Co = 0.5 in the simulations, corresponding to the overdamped regime. To 
determine the nonequilibrium voltage across the system, an external current / is injected 
uniformly with density J = I / L along one edge of a square array of size Lx L and extracted 
at the opposite one. Periodic boundary conditions are used in the transverse direction to 
the current. The average voltage drop V across the system is given by 

and the average electric field hj E = V/L. The dynamical equations were integrated nu- 
merically using typically a time step 6t = 0.03r (r = Ti/ {2eRoIc)) and averages computed 
with 2 X 10^ time steps for each calculation. Lattices of sizes L = 21 and L = 34 were used 
in the simulations without significant size dependence. 

In Fig. 1, we show the nonlinear I-V characteristics as function of temperature. At 
small current densities J, the nonlinear resistivity E/J shows a linear contribution where 
E/J is a constant. The linear resistivity decreases rapidly with decreasing temperature. For 
increasing J, the resistivity cross over to a nonlinear behavior at a critical current which also 
decreases with decreasing temperature. For T < 0.25, the linear behavior presumably occurs 
at current densities much smaller than the lowest value used in the simulations, J = 0.02. 
However, in order to confirm this behavior a much better statistics and long equilibration 



times would be required to obtain reliable data. This is consistent with a kind of frozen 
state for temperatures below T ~ 0.25 as indicated in Monte Carlo simulations! where a 
nonzero Edwards-Anderson order parameter q{t) =< Si where S = (cos(^^), sin(^^)), 
averaged over long times t, was observed below this temperature. This implies a diverging 
relaxation time r ~ /q°° q(t)dt which prevents one to obtain fully equilibrated data in this 
low temperature regime. 

In Fig. 2 we show the behavior of the linear resistance as a function of temperature, 
estimated from E/ J ai the smallest current. It appears consistent with an activated Arrhe- 
nius behavior with an estimated energy barrier Ei, ~ 0.45. We note that this energy barrier 
is higher than the single vortex barrier in a periodic array0, Ei, ~ 0.19. This suggests that 
many vortices are involved in the activation process and the vortex dynamics could be the 
result of collective motion of correlated vortices within a finite length scale. 

The behavior of the nonlinear I-V characteristics described above can be understood 
within a scaling theory which assumes a second-order phase transition at zero temperature^. 
This theory has been applied to the study of models of vortex glasses in disordered two- 
dimensional superconductorji0 but it can be equally applied to our case once a zero tem- 
perature transition is accepted. Basically, the behavior of the I-V curves is strongly affected 
by the increasing correlation length ^ which is assumed to diverge at zero temperature as 
^ T^" with a thermal correlation exponent u. At small current densities J, the behavior 
is dominated by vortex dynamics at length scales larger than ^ where the system behaves 
as a vortex liquid with a linear resistance Rl. At higher currents, one probes length scales 
smaller than ^ where the vortex dynamics is expected to be nonlinear. Near the transition, 
the nonlinear I-V characteristics can be cast into a scaling form@'i 

^ 9(7^) (4) 



JRl 

where g is a. scaling function which has the property that g{x) —>■ 1 when x 0. Since at 
nonzero temperatures the correlation length is finite, vortex motion will proceed by thermal 
activation leading to a linear resistance in the form 



Rl oc exp{-Eb/T) (5) 

where Eh is a typical barrier height for vortex motion within a correlation length ^. In 
general, the barrier energy should also scale with ^ as -Efe ~ ^"^ with a new exponent ^Z' > 
but the present data do not allow a clear identification of this behavior. The Arrhenius plot 
in Fig. 2 is consistent with the simplest scenario ip = ^ with Ef, ~ 0.45. The characteristic 
time scale r should also have an exponential behavior r ~ exp{Eh/T) leading to very 
slow relaxation at low temperatures consistent with the results above and Monte Carlo 
simulation^. The estimated Ef, being larger than the single vortex pinning potentiallil 
Ep = 0.19 is also consistent with a relaxation process involving collective rearrangement 
of pinned vortices. From the scaling behavior above, one also sees that the characteristic 
current density at which nonlinear behavior is expected to set in varies as Jni ~ T^^'^. In 
a finite system, there will also be a dependence on the system size but the system sizes we 
have used, L = 21 and L = 34, do not show significant corrections. 

In Fig. 3, we show a scaling plot of the I-V characteristics data of Fig. 1. From the 
above scaling analysis all data should collapse on to the same curve if u is chosen correctly, 
assuming size corrections are small. Choosing different values of u, we find that a reasonable 
scaling behavior is obtained for u ~ 0.9 ± 0.2. This value of u should be compared with the 
corresponding value for the gauge glass modeli of disordered superconductors in a magnetic 
field, z/ ~ 2. As expected, the different nature of the two vortex glass transitions is reflected in 
the values of this critical exponent. The sharp difference between the critical exponents also 
leads to a completely different behavior for the crossover current density which is expected 
to behave as Jni ~ in contrast with the ~ behavior for models of vortex glasses in 
disordered superconductors. In the experimental systems additional complications may arise 
since a higher order rational value can always be found within the confidence interval of a 
measured value of /. 

In conclusion, the current- volt age characteristics of an ordered Josephson junction array 
at an irrational fiux quantum per plaquette was studied using an scaling analysis of a zero 
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temperature transition. The correlation critical exponent is found to be significantly dif- 
ferent from the corresponding value for vortex-glass models in disordered two-dimensional 
superconductors. This leads to a current scale where nonlinearities sets in decreasing with 
temperature roughly as in contrast with the behavior expected for disordered models. 
It is suggested that this behavior could provide an experimental signature for an array at 
irrational frustration. 
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FIGURES 

FIG. 1. Nonlinear current- voltage characteristics at different temperatures for a system size 
L = 34. 

FIG. 2. Temperature dependence of the linear resistance as estimated from E/J at J = 0.02. 
FIG. 3. Scaling plot of the nonlinear current-voltage characteristics of Fig. 1 for z/ = 0.9. 



9 




Figure 1 
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Figure 3 
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